The effect of large-decoherence on 
mixing-time in Continuous-time quantum 
walks on long-range interacting cycles 

a^ : 

O ■ S. Salimi0, R. Radgohar 1 

. Faculty of Science, Department of Physics, University of Kurdistan, 

Pasdaran Ave., Sanandaj, Iran 

O 

o 

o 

m 

> ■ 

3 ' Abstract 

O^, In this paper, we consider decoherence in continuous-time quantum walks 

on long-range interacting cycles (LRICs), which are the extensions of the 
cycle graphs. For this purpose, we use Gurvitz's model and assume that 
1^ , every node is monitored by the corresponding point contact induced the 

Q"^ I decoherence process. Then, we focus on large rates of decoherence and 

' calculate the probability distribution analytically and obtain the lower and 

upper bounds of the mixing time. Our results prove that the mixing time 
is proportional to the rate of decoherence and the inverse of the distance 
Q>^ ' parameter (m) squared. This shows that the mixing time decreases with 

CD ■ increasing the range of interaction. Also, what we obtain for m = is 

in agreement with Fedichkin, Solenov and Tamon's results [48] for cycle, 
K/i ' and see that the mixing time of CTQWs on cycle improves with adding 

}_j ■ interacting edges. 
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1 Introduction 



Random walks or Markov chains on graphs have broad apphcations in var- 
ious fields of mathematics |1J, computer sciences [2] and the natural sci- 
ences, such as modeling of crystals in solid state physics [3]. The quantum 
walk(QW) is a generalization of the classical random walk(CRW) developed 
by exploiting the aspects of quantum mechanics, such as superposition and 
interference [HIS]. The QW is widely studied in many distinct fields, such as 
polymer physics, solid state physics, biological physics and quantum com- 
putation [SIEIE]- There are two distinct variants of QWs: the discrete-time 
QWs(DTQWs) [a do] and the continuous-time QWs(CTQWs) [H]. In the 
CTQW, one can directly define the walk on the position space, whereas 
in the DTQW, it is necessary to introduce a quantum coin operation to 
define the direction in which the particle has to move. Experimental im- 
plementations of both QW variants have been reported e.g. on microwave 
cavities p^, ground state atoms ^13], the orbital angular momentum of pho- 
tons [11] , waveguide arrays [15] or Rydberg atoms [16l[T7] in optical lattices. 
Recently, quantum walks have been studied in numerous publications. For 
instance, the DTQWs have been considered on random environments |18j . 
on quotient graphs [19], in phase space [20] and for single and entangled 
particles [21] • Also, the CTQWs have been investigated on the n-cube [22], 
on star graphs [231 |2l], oii quotient graphs [25], on circulant Bunkbeds [26] . 
on trees [27], on ultrametric spaces [28], on odd graphs [29] , on simple one- 
dimensional lattices [30] . via modules of Bose-Mesner and Terwilliger alge- 
bras [31], via spectral distribution associated with adjacency matrix [32] and 
by using krylov subspace-Lanczos algorithm |33| . 

In this paper, we focus on the CTQWs on the one-dimensional networks. 
These networks provide the better understanding of the physical systems 
dynamics. For examples, they have been used to describe the behavior of 
metals in solid state physics [3] , to explain the dynamics of atoms in optical 
lattices \13\ [T6] , to demonstrate Anderson localization in the systems with 
energetic disorder ^34j and to address various aspects of normal and anoma- 
lous diffusion [35] . The simplest structure of the one-dimensional networks 
generated by connecting the nearest neighbors can be modeled the physical 
systems with short-range interactions. Gravity and electromagnetism are 
the only known fundamental interactions extending to a macroscopic dis- 
tance. Due to its basic importance, it has been along tradition to search 
extra long-range interactions [36] . In some of the above mentioned experi- 
ments, e.g. the clouds of ultra-cold Rydberg atoms assembled in a chain over 
which an exciton migrates, the trapping of the exciton occurs at the ends of 



2 



the chain, one finds long-range interactions have to be considered [n\ I37j . 
Recently, it has been shown, from the view point of Quantum Electrodynam- 
ics, the CTQW for all long-range interactions in a linear system decaying 
as R~'^ (where R is the distance between two nodes of the network) belong 
to the same universality class for ly > 2, while for classical continuous-time 
random walk (CTRW) universality only holds for u > 3 [38j. It has been 
shown the long-range interaction lead to a slowing-down of the decay of the 
average survival probability, which is counter intuitive since for the corre- 
sponding classical process one observes a speed-up of the decay [39]. In [30], 
authors studied coherent exciton transport on a new network model, namely 
long-range interacting cycles. LRICs are constructed by connecting all the 
two nodes of distance m on the cycle graphs. Since all the LRICs have the 
same value of connectivity k = 4 (the number of edges which exit from ev- 
ery node), thus LRICs provide a good facility to study the influence of long 
range interaction on the transport dynamics on various coupled dynami- 
cal systems, including Josephson junction arrays [41j, synchronization |42j . 
small-world networks [43] and many other self-organizing systems. All of the 
mentioned articles have focused on a closed quantum system without any 
interaction with its environment. Recently, more realistic analysis of quan- 
tum walks by using decoherence concept has been provided on line |44j . 
on circulant [l5l [l6] and on hypercube [U]. In [JH], the effect of decoher- 
ence in the CTQWs on cycles has been considered analytically. The results 
showed that, for small rates of decoherence, the mixing time improves lin- 
early with decoherence, whereas for large rates of decoherence, the mixing 
time deteriorates linearly towards the classical limit. In [l9], the authors 
have studied the influence of small decoherence in the CTQWs on long- 
range interacting cycles (LRICs). The authors proved that the mixing time 
is inversely proportional to the decoherence rate and also independent of the 
distance parameter m. Moreover, they showed the mixing time upper bound 
for CTQWs on LRIC, remain close to its value in the absence of shortcut 
links(cycle). Now, we want to investigate the effect of large decoherence 
on the mixing time in the CTQWs on LRICs. For this end, we use of an 
analytical model developed by Gurvitz [50]. In this model, every vertex is 
monitored by an individual point contact induced the decoherence process. 
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Figure 1: Long-range interaction cycles G(8,3) and G(10,4) 



We calculate the probability distribution analytically then, for large rates 
of decoherence, obtain the lower and upper bounds of the mixing time. Our 
analytical results prove the bounds of the mixing time are proportional to 
the decoherence rate that in agreement with [48j's results . Moreover, we 
prove that these bounds are inversely proportional to square of the distance 
parameter m, i.e. the mixing time decreases with increasing m. 
This paper is organized as follows: In Sec. 2, we briefly describe the network 
structure LRIC. The CTQWs on LRICs are considered in Sec. 3. In Sec. 4, 
we study the effect of decoherence in CTQWs on LRICs, then we focus on 
the large rates of decoherence and calculate the probability distribution in 
Sec. 5. In Sec. 6, we define the mixing time and obtain its lower and upper 
bounds. We conclude with a summery of results in Sec. 7. 

2 Structure of LRIC 

Long-range interacting cycle (LRIC) can be generated as follows |40j : First, 
the network be composed of a cycle graph. Second, two nodes of distance 
m on the cycle graph are linked by an additional bond. We continue the 
second step until all the two nodes of distance m have been connected. LRIC 
denoted by G{N, m), is characterized by the network size N and the long- 
range interaction parameter m. Fig. 1 shows the sketches of G{8, 3) and 
G(10,4). 

3 CTQWs on LRICs 

In general, every network is characterized by a graph whose bonds connect 
nodes with a wide distribution of mutual distances. Algebraically, every 
graph corresponds with a discrete Laplace operator A. We introduce the 
states \j) which are localized at the nodes j of the graph and take the set 
{\j)} to be orthonormal. We assume that transition rates 7 between all 
nodes are equal and 7 = 1. The non-diagonal element of matrix A (Aij) 
equals 1 if nodes i and j are connected by a bond and otherwise. The 
diagonal element An equals —ki that ki is the number of bonds which exit 
from node i 151] . Since the states |j) span the whole Hilbert space, 
the time evolution of a state \j) starting at time is determined by the 
system Hamiltonian H = A as \j,t) = U{t)\j), where U{t) = exp[—iHt] 
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is the quantum mechanical time evolution operator |1H, I51| . Thus the 
Hamiltonian matrix H of G{N, m) (m > 2) can be written as 



H, 



1-] 



(i\H\j) 



-4, ifi = i; 

1, \ii=j± 1; 

1, if i = J lb m; 

0, Otherwise. 



The Hamiltonian acting on the state \ j) is given by 

H\3) = -4|i) + \j - 1) + |j + 1) + \j -m) + \j + m). 



(2) 



The above equation is the discrete version of the Hamiltonian for a free 
particle moving on the cycle. Using Bloch function [3l [52] in solid state 
physics, the time independent Schrodinger equation can be written as 

H\ll,n)=En\ll>n). (3) 

The Bloch states \ipn) can be expanded as a linear combination of states |j) 



1 ^ 



(4) 



Substituting Eqs. (2) and (4) into Eq. (3), we obtain the eigenvalues of the 
system as 



Er, 



-4 + 2cos(6'„) + 2cos(m6'„). 



(5) 



The Bloch relation can be obtained by projecting \'4)n) on the state |j) such 
that Vn(j) = {Mn) = e-'^^"^^/N. It follows that On = 2m: /N with n 
integer and n G [0, A^). From Schrodinger equation, we have 

if^j^\Mt)) = H\Mt))- (6) 

By assuming h = 1 and |^n(0)) = |0), the solution of the above equation is 
\Mt)) = e-'"'\0). 

The probability to find the walker at node j at time t is given by 
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4 The Decoherent CTQWs on LRICs 



Here, we want to study the appearance of decoherence in the CTQWs and 
obtain the probabihty distribution P{t). For this end, we make use of an 
analytical model developed by Gurvitz \50\ I53j . In this model, a ballistic 
point-contact is placed near each node of network that is taken as non- 
invasive detector. Gurvitz demonstrated that measurement process is fully 
described by the Bloch-type equations applied to whole system. These equa- 
tions led to the collapse of the density-matrix into the statistical mixture in 
the course of the measurement process. 

According to Gurvitz model, the time evolution of density matrix for our 
network (LRIC), has the following form 

^Pj,fc(0 = ■l[—Pj-l,k — Pj+l,k — Pj-m,k — Pj+m,k + Pj.k-1 

(7) 

+Pj,k+1 + Pj,k-m + Pj,k+m] — r(l — dj^k)Pj,k, 

that density matrix p{t) is as p{t) = \^p{t)){■^p{t)\ and T is the decoherence 
parameter. Also, we observe that Pjit) = pjj{t). 
We define the variable Sj^k as [48] 

Sj,k = i^~^Pj,k (8) 
and by substituting it into Eq. (7), we obtain 

do 1 r C _i_ Q m+1 c „'m+l q c 

dt — 41- '-'j — l,k ~r i-^j+l,k ' ^j—m,k ' »-'j+m,fc '^j,k—l 

(9) 

+ Sj^k-m + i ^ S'j^fc+m] — r(l — 5j^k)Sj^k- 
Note that by assuming m = 0, we achieve the relations mentioned in |48j . 

5 Large Decoherence 

Firstly, we review the results obtained in [49j for the decoherent CTQWs 
on LRICs with small rate of decoherence. The mixing time upper bound for 
the odd values of m is as Tmix{() < ^ ^"^^ even ms is as 

Tmix{^) < r ^'^(t')[i^t]' Thus, the mixing time upper bound for odd m is 
larger than the mixing time upper bound for even m. In addition to these 
relations show the upper bound of the mixing time is inversely proportional 
to decoherence rate F and independent of the distance parameter m. 
In the following, we assume that the rate decoherence F is large (F S> 1) 
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and consider its effect in CTQWs on LRICs. 

Firstly, as mentioned in [l8], we define the diagonal sum as 



AT-l 

Dk = ^ Sj j_^_i.(^jyiodN)^ (10) 
j=0 

where the indices are treated as integers modulo A'^. 
From Eq. (9), one can achieve the following form 

j^Dk = -r{l-6k,o)Dk. (11) 

Thus the minor diagonal sums {D^ for k ^ 0) decay with characteristic time 
of order l/F. 

Also, by applying Eq. (9) for the elements on the two minor diagonals 
nearest to the major diagonal, we observe that these elements make limit to 
small values of order l/F. 

Using the above way for the other diagonals, it is yielded the secondary set 
of matrix elements in terms of nearness to major diagonal will be of the 
order of 1/F^, etc. 

Now, we consider only matrix elements of order of l/F and get 



— Q. . -, _|_ Q. . , , _|_ Q. . -i_A-m+lq, 1 



(12) 



■^iJ+i ~ ll'S'j+ij+i - Sjj] - rSjj+i, (13) 



•^ij-i ~ ji~^3-^,3-i + Sjj] - rSjj-i, (14) 



^j,3+-m — ^[ '^"^^'^^j+m,j+m + i"^^^Sjj] TSjj+m (15) 



^3,j-m ~ ^[ ^ ^j-m,j-m + i Sjj] TSjj-m- (16) 
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For simplicity's sake, we define 

After some algebra, one can get 

<^'j = \l-dj-i + dj+fj-fj-m], (18) 

dj = ^[aj+i-aj]-'^dj, (19) 



/j = i[a,+„-a,]-r/,-. (20) 
Differentiation of the above equation gives 

< = ^[-4_i + 4 + /j-/j_J, (21) 



< = ^K.+i-a;.]-r4, (22) 

f!=l^+m-^j]-'^f^- (23) 

Let us now assume that we have the solutions of Eqs. (21), (22) and 
(23) as 



N-l N-1 N-1 

aj = J2 A^e^er^^\ dj = ^ Dke'-^e-^>'\ fj = i^fee^e-^'^*, (24) 

fe=0 k=0 k=0 

that 7fc, Ak, Dk and Fk are unknown. 
To obtain 7;^, we substitute these solutions into the above mentioned equa- 
tions. 

All + i?(l(-eT + 1)) + F(^(-eT^ + 1)) = 0, 

A(^(e^-l))+i5(7|-7fcr) = 0, (25) 
A(^(e^-l)) + F(72-7fer) = 0. 
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Note that for having nontrivial solutions, the determinant of the coeffi- 
cients matrix must be zero. Thus, we get 



llill - lkT)bl - 7.r + ^(sin2(^) + sin^(^))] = 0, 



AT 



and achieve 



7fe = < 



7M=2r(sin'(f)+sm2(^))> 
7M=r-2^(sin2(^) + sin2(l^)), 

7fc,2 = 0, 

I 7M = r- 



(26) 



(27) 



The general solutions of Eq. (21), (22) and (23) are as 

N-l 



fe=o 



N-l 



fe=o 



N-l 



fe=o 

2iTi 



where lo = e~T^ . We replace 7^3 in Eqs. (25) and use of the initial conditions 
Qj = 5jfl and dj = for j = 0, A/" — 1. We set the constant coefficients 
into the others and obtain 



Afc,o^l, Afe,i:^2j4(sin2(^)+sin2(^)), ^^,3 = 0, 



(31) 



Dk,o ^ f sin(f ) exp(^), Dk,i ^ ^ sin(f ) exp(^). 



Trfe ^ 



(32) 



F,,o ^ f sin(^) exp(Hf-), ^ sin(^) exp(Hf-), 



f iirkm 



km ^ 



rkm ^ 



i nkm iirkm 
Fk,3 - ^ sm(-^) exp(-^). 



(33) 
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The full solution for S{t) has the form as 



if j=k; 
dj/2, iij = k±l; 
fj/'^, iij = k±m; 
0, otherwise. 



(34) 



Thus, the probability distribution Pj{t) is the same aj. At large rates of 
decoherence F, Eq. (28) reduce to 



N-l 



k=0 



irkm 
N 



(35) 



One can see the above distribution do not converge to any stationary 
distribution. The reason being that the evolution of the QW, as mentioned 
in Sec. 3, is given by the unitary operator U = e~^^^ . This is the fact 
that unitary operators preserve the norm of states, and hence the distance 
between the states describing the system at subsequent times does not con- 
verge to zero[9]. This implies that the probability distribution of CTQW 
does not converge to CRW. 



6 The bounds of mixing time 

Mixing time is the time it takes for the walk to approximate the uniform 
distribution [Ml [9], i.e. 

N-l ^ 

Tmix = min{T : ^ |P^.(t) - - | < e}. (36) 
i=o 

From Eq. (35), we obtain 

N-l N-l N-l 111 1 

-n\ = Y.\nY. -M^isin\-) + sin\^))t]e^ - -|,(37) 

j=0 j=0 k=0 

that simplifies to 

E W - ^1 = ^ E I E exp[^(sin^(^) + sin^(^))t] cos(^)|.(38) 

1=0 7=0 k=l 
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Lower bound 

To obtain a lower bound, we apply the following way: 
We retain only the terms j = and k = 1, N — 1. 

V^'l m 1|> I m i| 1 V^' , -(sin^(^)+sin^(^))t . 
}_^\aj{t)--\> |ao(t)-^| = ^ 2^exp( L_ ^ 

j=0 fc=l 

2 (sin^(^)+sin^(2^))t ^ 

- iv""P^ 2r 

According to the mixing time definition, we have 

2T 2 



(39) 



and for iV > 1 



2riV^ 2 
7r^(l + m^j iVe 



Note that for m = 0, 



2TN'^ 2 

Tzo^,, ~ ^^ln(— ), (42) 



that is the same (48j's result. One observe the mixing time lower bound in 
CTQWs on cycle decreases with adding interacting links. 

Upper bound 
Now, we want to obtain upper bound as following 

E = ^ E I E ^M^{s^n\ — ) + sin2(^))t] cos(^: 

i=o j=0 k=l 

< ^ E E exp( (-^M^)+sm(^)) ^ 

j=0 A:=l 
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(44) 



Since for < x < ^, there is sinx > [5l], we get 

Af-l ^ N~l[N/2] _i ^ 

5^ |a,(t) - -I < f ^ E exp(— (sin2(^) + sin\^))t) 

j=0 j=0 k=l 

< # E E «^p(i^( — ^2 — )t) 

j=0 k=l 

9 'v^' —1 ,2k + 2km. . 

< I E E ^^P(ir( )0 

9 ^^^^^ , — 1 ,2A; + 2A;m, , 

< fEE^^P(ir( ^2 )^) 

that in the third inequahty, we used of the relation k'^ > k for k > 1. 
^-1 1 2 



Thus, 

riV^ ^ ^2 + e 
2(1 +m2 

Note that for m = 0, we get 



Tupper = , ln( — — ). (46) 



riV2 2 + e 
Tapper = ^^ln( ). (47) 

which is in agreement with [38] 's result. One can see that the mixing time 
upper bound for LRIC is smaller than the mixing time upper bound for cycle. 
Based on the above analysis, the mixing time in CTQWs on cycle decreases 
with adding interacting links. Moreover, we observe that the mixing time is 
proportional to the square of distance parameter m such that it decreases 
with increasing the interaction range. Now, we want to compare the lower 
bound with the upper bound. Since Eq. (41) was obtained by assuming 
^ 1 and e <C 1, we have In(-j^) < In(^i^). Thus, we can write 

2rA^2 2 , ^ riV2 2 + e, 

■ ln( — ) < r„i, < — ln( ). (48) 



^2(1 + m2) 'Ne' ™" 2(1 +m2) 
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7 Conclusion 



We studied the continuous-time quantum walk on long-range interacting 
cycles (LRICs) under large decoherence F » 1 . We obtained the probability 
distribution analytically and found the mixing time is bounded as 

■ ln( — ) < T„i, < -— In(^). (49) 



^2(1 + m2) ^^Ne' ™" 2(1 +m2) 

We proved that the Tmix is proportional to decoherence rate F that 
this result accord with the conclusion is obtained in [38]. In [l9] has been 
proved that Tmix-, for small rates of decoherence, is independent of distance 
parameter m, while in this paper we showed that Tmix, for large rates of 
decoherence, is inversely proportional to square of m. In other words, the 
mixing time decreases with increasing the range of interaction. Also, since 
LRIC is the same generalized cycle, by replacing m = in all the above 
relations, one achieves p[8]'s results. Moreover, we proved the mixing times 
improve with adding interacting links. 
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